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These lecture notes focus on the application of ideas of locality, in particular Lieb-Robinson 
bounds, to quantum many-body systems. We consider applications including correlation decay, 
topological order, a higher dimensional Lieb-Schultz-Mattis theorem, and a nonrelativistic Goldstone 
theorem. The emphasis is on trying to show the ideas behind the calculations. As a result, the 
proofs are only sketched with an emphasis on the intuitive ideas behind them, and in some cases we 
use techniques that give very slightly weaker bounds for simplicity. 

This is a preliminary version of the lecture notes, with the goal of getting the notes out close to 
' the end of the school. Comments welcome. 
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^ ■ I. INTRODUCTION AND NOTATION 

. . . 

^ The basic problem studied in quantum many-body theory is to find the properties of the ground state of a given 

Hamiltonian. Expressed mathematically, the Hamiltonian H is a, Hermitian matrix, and the ground state, which we 
write "00, is the eigenvector of this matrix with the lowest eigenvalue. The properties we are interested in studying 
are expectation values of various observables: given a Hermitian matrix O, we would like to compute the expectation 
value {ipo,Oipa)- A note on notation: we write an inner product of two vectors as {v,w), rather than the notation 
I— I {v\w), as the first notation is more common in math and also will be more clear given the large number of absolute 
i-S^ , value signs we will also use later). Sometimes we will write the expectation value in the ground state simply as (O). 

The above paragraph sets out some of the mathematics of quantum mechanics, but it must seem quite dry. Someone 
, reading that who is not familiar with quantum theory could be excused for thinking that the problem is essentially one 
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of linear algebra, and that the basic tool employed by quantum physicists is a linear algebra package to diagonalize large 
matrices. In fact, while the method of "exact diagonalization" on a computer is an important technique in studying 
^ ^ ^ ■ quantum systems, it is only a small part of how physical problems are studied. The feature that distinguishes the study 
of quantum many-body systems is locality of interactions. Consider a typical Hamiltonian, such as the one-dimensional 
,—1 . transverse field Ising model: 
> • 

■ N-l N 

m: h = -jJ2s!sUi + bY,s^- (1) 



in 



1=1 



The very notation we employ to describe this Hamiltonian implicitly assumes local interactions. To be more precise, 
throughout these notes, we have in mind quantum systems on a finite size lattice. We associate a D dimensional Hilbert 
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• space with each lattice site, and the Hilbert space of the whole system is the tensor product of these spaces. We use N 
[ to represent the number of lattice sites, so that the Hilbert space on which a Hamiltonian such as (U) is defined is 
'' • dimensional. A term such as Sf,^]^ is a short-hand notation for the term /i (g)/2(8) ...®/i-i (^S'f (SJS'f^.]^ ®/i+2 (8)... (8)/iv, 
' where Ij is the identity operator on site j. This local structure of the interaction terms will greatly constrain the 
. properties of ipo, in particular if there is a spectral gap. 
^ ' We will usually use symbols i,j,k,... to denote lattice sites and we use letters X,Y,Z,... to denote sets of lattice 
sites. We use A to denote the set of all lattice sites. We use \X\ to denote the cardinality of a set X. 

We say that an operator O is "supported on set A" if we can write O as a tensor product of two operators: 

O = Ia\a ® P, (2) 

where Ia\a is the identity operator on the sites not in set A (note that A \ ^ denotes the set of sites not in A) and 
P is some operator defined on the Z?'"*' dimensional Hilbert space on set A. For example, the operator SfSfj^i is 
supported on the set of sites {i, i + 1}. CoUoquiually, instead of saying that O is supported on set A, one often hears 
that O "acts on set ^" , although we will avoid that terminology. 

Finally, we define ||0|| to represent the "operator norm" of an operator O. If O is Hermitian, the operator norm is 
equal to the absolute value of the largest eigenvalue of O. For arbitrary operators O, we define 

||0|| = max^j.0|=i 10-01, (3) 

where the maximum is taken over vectors -0 with norm 1, and \Oili\ denotes the norm of the vector O0. 

Using two simple assumptions, that the Hamiltonian has local interactions and that the Hamiltonian has a spectral 
gap, we will be able to prove a wide variety of results about the ground state of the Hamiltonian. We will consider two 
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different cases of a spectral gap. In one case, we consider a Hamiltonian with a unique ground state, ipo, which has 
energy Eq, with the next lowest energy state having energy Ei with Ex — Eq> AE. Here, AE is called the "spectral 
gap" . In the other case, we consider a Hamiltonian with several degenerate or approximately degenerate low energy 
states, ipQ, for a = l,...,k. Here, k is the number of different low energy states. Let these states have energies Eq. 
Then, assume that there is a gap AE separating these ground states from higher energy states. That is, the k + 1-st 
smallest eigenvalue of the Hamiltonian is greater than or equal to Eq + AE for all a = 1, fc. 

We will begin by considering properties of correlations in these systems, and prove an exponential decay of corre- 
lation functions. We then consider how the ground state of the system changes under a change in the Hamiltonian, 
and use this to prove a non-relativistic variant of Goldstone's theorem. Finally, we apply these techniques to more 
interesting systems with topological order, beginning with the Lieb-Schultz-Mattis theorem. However, before we can 
do any of this, we need to more precisely define the locality properties of the Hamiltonian and to prove a set of bounds 
on the propagation of information through the system called "Lieb- Robinson bounds" . 

II. LOCALITY AND LIEB-ROBINSON BOUNDS 

To define locality, we need a notion of distance. We introduce a metric on the lattice, dist(i, j). For example, on a 
square lattice one may consider a Manhattan metric or a Euclidean metric. In one dimension on a system with open 
boundary conditions, it is natural to consider the metric dist(i, j) = \i — j|, while for a one dimensional system with 
periodic boundary conditions and N sites it is natural to consider the metric dist(i, j) — min„|i — j + nN\, where the 
minimum ranges over integers n. We define the distance between two sets A, B to be 

dist(A, B) = miuig^ jgBdist(i, j). (4) 

We define the diameter of a set A by 

diam(y4) = maxi jgAdist(i, j'). (5) 

We consider Hamiltonians 

H = Y^Hz, (6) 
z 

where Hz is supported on a set Z . We will be interested in studying problems in which decays rapidly with the 

diameter of the set Z . For example, in the one-dimensional transverse field Ising model Hamiltonian ([T]), with metric 
dist(i, j) = |« — il, the only terms that appears have diameter or 1 (the magnetic field term has diameter zero and 
the Ising interaction has diameter one), so \\Hz\\ = for diam(Z) > 1. These interactions are "finite range"; we will 
also be able to consider cases in which the interactions have a slower decay, such as an exponential decay (one can 
prove various results about power law decay also, which we will not consider in these notes). 

Note that the reason one chooses one particular metric over another is to make the terms in the Hamiltonian local 
with respect to that metric. Note that Hamiltonian ([T]) has open boundary conditions, in that site iV does not interact 
with site 1. Instead, if we considered a one-dimensional transverse field Ising 

Af-l N 

H = - jy^^ SjS^^i - JS%Sl + By^Sj, (7) 

1=1 i=l 

with periodic boundary conditions, the term S^Sf would have diameter A'^ — 1 with respect to the metric dist(i,j) — 
\i — j\. However, if we instead consider the metric dist(i, j) — min„|i — j + nN\ suggested above, then all of the 
interactions have diameter or 1 again. 

As an aside, one might wonder why we do not just pick the metric dist(i, j') — for all In this case, all of the 
interactions in the Hamiltonian have diameter zero and are therefore "local" . However, the kinds of results we will 
prove below have to do with exponential decay of various quantities with distance (for example, decay of correlations 
as a function of spacing between operators), and all of the results would become trivial in the case of the metric 
dist(i, j) = 0. Thus, we want to pick a metric such that our Hamiltonian is local, but such that the set of sites A has 
a large diameter. 

It is not necessary that the system be defined on a regular lattice. For example, we could imagine a graphical 
structure to the interactions. Suppose each site is associated with the vertex of a graph, and the Hamiltonian is a 
sum of terms Hz where each Z contains only two sites and Hz is non- vanishing only if there is an edge of the graph 
connecting the sites in Z. Then, the shortest path metric on the graph gives us a metric for which Hz is non-vanishing 
only if Z has diameter or 1. 
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We now consider the Lieb-Robinson bounds. These are bounds describing time evolution of operators in a local 
Hamiltonian. While most of our interest in these notes is in static properties of systems (such as correlation functions 
at a given instant in time), it turns out that an understanding of the dynamical properties as a function of time is very 
useful to prove results about the statics. The bounds were first proven in Then, a proof which does not involve 
the dimension of the Hilbcrt space dimension on a given site was given in Q . More general lattices were considered in 
but only using a proof which depended on the Hilbert space dimension again. The dimension independent proof 
for arbitrary lattices was given in '3l| and this is what we follow. 

Various sorts of Lieb-Robinson bounds can be proven. We consider first the case of exponentially decaying interac- 
tions in some detail, then we consider either kinds of decay in the interaction (either finite range interactions or other 
types of decay). We define the time-dependence of operators by the Heisenberg evolution: 

0{t) = exp{iHt)Oexp{~iHt). (8) 

Then, 

Theorem 1. Suppose for all sites i, the following holds: 

^ ||i/jf|||X|exp[/idiam(X)] < s < oo, (9) 

for some positive constants /i, s. Let Ax, By he operators supported on sets X, Y , respectively. Then, i/dist(X, Y) > 0, 

\\[Axit),BY]\\ < 2Px||||By||^exphAidist(i,y)] [e2^l*l-l] (10) 

< 2\\Ax\\\\BY\\\X\exp[~fidi8t{X,Y)] [e^-l*! - 1 

Before proving the theorem, let us describe the physical meaning of this theorem. Suppose Ax is some operator. 
Let Bi{X) denote the ball of radius I about set X. That is, Bi{X) is the set of sites i, such that dist{i,X) < I. 
Following define 

Axit) = J dUUAxm\ (11) 

where the integral is over unitaries supported on the set of sites A \ Bi{X) with the Haar measure. Then, is 
supported on Bi{X). Since UAx{t)W = Axit) + U[Axit), W], we have 

\\A^xit)-Ax{t)< J dU\\[Ax,U]\\. (12) 

Using the Lieb-Robinson bound ([TUl) to bound the right-hand side of the above equation, we see that (t) is expo- 
nentially close to Ax if I is sufficiently large compared to 2si//i. Thus, to exponential accuracy, we can approximate 
a time-evolved operator such as Ax{t) by an operator supported on the set Bi{X). That is, the "leakage" of the 
operator outside the light-cone is small. 

Another comment regarding the assumptions on the theorem: because the term exp[/idiam(Z)] grows exponentially 
in Z, we need the norm of the terms Hz to decay exponentially in Z. On finite dimensional lattices (such as a 
hypercubic lattice), the cardinality of Z, \Z\, is bounded by a power in the diameter of Z. So, on such lattices, if 
the terms Hz in the Hamiltonian decay exponentially in diam(.Z), we can find a fi such that the assumptions of the 
theorem are satisfied. By expressing the theorem as we have, however, the theorem can be applied to models defined 
on, for example, arbitrary graphs, where the cardinality of a set might not be bounded by a power of its diameter. 

One final remark: we have implicitly assumed that all of the operators are bosonic, in that operators supported 
on disjoint sets commute with each other. One can straightforwardly generalize all of this to the case of fermionic 
operators also, so that two fermionic operators which are supported on disjoint sets anti-comniutc but that two bosonic 
operators, or one bosonic and one fermionic operator, commute on disjoint sets. In this case, one can instead prove a 
bound on the anti-commutator of two fermionic operators at different times. 



A. Proof of Lieb-Robinson Bound 



The proof of the Lieb-Robinson bound we now give can be straightforwardly adapted to time-dependent Hamilto- 
nians, though we only present the proof in the time-independent case for simplicity of notation. 
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Recall that we assume that |A| is finite. If it is necessary to consider the infinite volume limit, we take the limit 
after deriving the desired Lieb-Robinson bounds which hold uniformly in the size of the lattice. Our essential tool is 
the series expansion (|27|) below for the commutator [A{t),B]. Let A be supported on X and B be supported on Y. 
We assume i > because negative t can be treated in the same way. Let e — t/N with a large positive integer TV, 
and let 

tn = for n = 0,l,...,N. (13) 

Then we have 

II [Ait),B]\l - II [A(0), B]\lJfe. \\mn,.).B]\\-\mt.),B]\l_ ^^^^ 

i=0 

In order to obtain the bound ([2^ below, we want to estimate the summand in the right-hand side. To begin with, 
we note that the identity, ||?7*0[/|| — \\0\\, holds for any observable O and for any unitary operator U. Using this 
fact, we have 

\\[A{t^+,),B]\\ - \\[Ait^),B]\\ = ||[A(e),i?(-i„)]|| - \\[AB(-t,M 

< \\[A + te[HA,A],B{-t,M - \\[A,B{-tn)]\\+0{e^) 
= \\[A + te[Ix,AlB{-tn)]\\ - II [A, S(-<„)]|| + ©(e^) 



(15) 



with 



z-.znx^di 

where we have used 

Aie) = A + ie[H^,A]+Oie^) (17) 

and a triangle inequality. Further, by using 

A + ie[Ix,A] = e^^-f^Ae-"^^ +0(e2), (18) 

we have 

\\[A + ie[Ix,A],B{-t,M < \\[e^-i-Ae-'''-,B{-t„)]\\+0{e^) 

= ||[A,e-"^^B(-i„)e"^^]|| +0(6^) 

< \\[A,B{-U)-te[Ix,B{-t,M+0{e^) 

< \\[A,B{-t^)]\\+e\\[A,[Ix,Bi^tr.m\+Oie^). (19) 

Substituting this into the right-hand side in the last line of (1151) . we obtain 

\\[A{tr^+,),B]\\-\\[A{tr^),B]\\ < e\\[A,[Ix,Bi-tr.m+0{e') 

< 2e\\A\\\\[Ix{tr.),B]\\+0{e'). (20) 

Further, substituting this into the right-hand side of (|14l) and using (|16l) . we have 

\\[Ait),B]\\-\\[A{0),B]\\ < 2\\A\\Y,ex\\[Ix{t,.),B]\\+0{e) 

n=0 

< 2\\A\\ E^x|l[^^(*«)'^]ll+^(^)- (21) 

Z:ZnX^(/) n=0 

Since Hz{t) is a continuous function of the time t for a finite volume, the sum in the right-hand side converges to the 
integral in the limit e ^ {N f c») for any fixed finite lattice A. In consequence, we obtain 



|[A(t),i?]||-||[^(0),i3]|| <2P|| r ds\\[Hz{s),B]\ 



(22) 



z-.znx^iH'- 
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We define 



CMX,0:= »p MM, ,23, 

AeAx ll^ll 



where Ax is the set of observables supported on the set X. Then we have 



CB{X,t)<CB{X,G) + 2 V \\Hz\\ I dsCB{Z,s 



(24) 



from the above bound ([2^ . Assume dist(X, F) > 0. Then we have Cb{X, 0) = from the definition of Cb{X, t), and 
note that 

Cb(Z,0)<2|1B||, (25) 

for Z n y =^ and 

Cb{Z,0))^0 (26) 
otherwise. Using these facts and the above bound ([M]) iteratively, we obtain 

"|t| 



CBiX,t) < 2 W^zAl f dsiCBiZ,,si) 

< 2 V \\HzA\ / dsiCB{Zi,Q) 

ry r, r^^.ru JO 



Zi:Zinx^<h 

+ 22 \\HzA\ Y \\HZ,\\ dsi t ds2CB{Z2,S2) 

Zf.Zinx^tli z.2:Z2r\ZiTt$ ■'^ •'° 

< 2\\Bmt\) E ii^^iii 

Zi:Zinx#0,ZinF#0 
Zi:Zinx#0 Z2:Z2nZi7t0,Z2ny#0 

+ 2iiBiiM! ^ iiHzAi E ii^^^ii E ii^^3ii + --- 

Zi:Zinx=^9 Z2:Z2n2i#0 23:23022^0, Z3ny#0 

(27) 

We now bound each term in Eq. (|27l) . using the assumption on the decay of terms in the Hamiltonian The first 
term is bounded by 2(2|t|) J2iex exp(— /^dist(i, Y). The second term is bounded by 

2^EEE E 1- (28) 

lex ZiBijeZi Z23j,Z2nY=^$ 

Recall that dist(i,y) < dist(«,j) +dist(j,y). Thus, exp[— ^dist(i, F)] exp[/idist(i, j)] exp[/idist(j, F)] > 1. Thus, the 
second term is bounded by 

2l^exp[-dist(i,y)]E E E '^'^P[^di«t(*,j)] E exp[Aidist(j,F)] (29) 

iexZiBijeZi Z23j,Z2nY^iS 

< 2^^exp[-dist(i,y)]E E E ^^P(^^i^*(*'^')) E exp[/^diam(Z2)] 

ieX Zi3i jeZi Z29j,22nY#0 

< 2^^exp[-dist(i,y)]E E E ^^Pt^^i'^*^*' 

iexZiBijeZi 

< 2-^^^exp[-dist(i,y)] E E l^i|exp[/^diam(Zi)]s 

ieX Zi3i 

< exph/idist(i,y)]2^^^s2. 
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Proceeding in this fashion, we bound the n-th term in the series ([27]) by 2 exp[— /idist(z, F)](2s|t|)"/n!. Adding 
these together, we arrive at the bound ([TU)) . 

B. The Lieb-Robinson Velocity and Finite- Range Interactions 

The bound (|10p is not in quite the most convenient form for later use. The most convenient form of the theorem 
is that, supposing the bound ([9]) holds, then there is a constant vlr depending only on s,/i such that for t < 
dist{X,Y)/vLii, we have 

\\[Ax{t),By]\\ < ^!l^gil)\X\\\Ax\\\\By\\, (30) 

where / = dist(X, Y) and g{l) decays exponentially in /. One may choose, for example, vlh = 4s//i and Eq. (pO)) will 
follow from Eq. (|10l) . Similarly, the operator Ax (t) can be approximated by an operator (t) supported on the set 

of sites within distance I — VLRt of the set X up to an error bounded by "^^•^*^ g(/)|X|||Ax||. 

The reader should note that in almost all applications of the Lieb-Robinson bound, the error term on the right-hand 
side of Eq. ([50]) associated with the "leakage" outside the light-cone is negligible. For example, in the next section, 
we will prove exponential decay of correlation functions. The proof will involve summing up various error terms. The 
most important error terms will not arise from the Lieb-Robinson bound. In fact, for clarity, the reader may wish to 
pretend that such error terms are zero on the first reading of any of the proofs in the rest of the paper, and only later 
worry about how large these errors are. 

Other types of interaction decay can be consider similarly. Suppose the interaction is finite range. For example, 
suppose Hz — for diam(Z) > R for some interaction range R. Then, we can find constants /i, s such that the 
assumption ^ is satisfied. However, we can do even better in the case of finite range interaction. Suppose each term 
in the Hamiltonian only has support on a set of two sites, and suppose we have R — I for simplicity and suppose 
\\Hz\\ < J for some constant J. We again use the series ((Tf)) . However, we note that we have now bounded CB{X,t) 
by a sum over paths on the lattice starting at sites in X and ending at sites in Y, with each path weighted by 
{2J\t\y/l\, where / is the length of the path. Such a weighted sum over paths is a well-studied problem in statistical 
physics. In this case, one finds that g{l) decays faster than exponentially in I (roughly, g{l) is exp(— const * as one 
may verify). Other cases of finite range interaction can be handled similarly. 

C. Other Types of Decaying Interactions 

One may also consider other types of decaying interaction, slower than exponential. This subsection will be useful 
in considering quasi-adiabatic continuation later, but may be skipped on first reading. 

The important things to remember from the calculations in the rest of the section are that one can treat decay 
other than exponential, and that decays which are slower than exponential decay give rise to Lieb-Robinson bounds 
with error terms that also decay slower than exponential. That is, we find error terms on the right-hand side of 
the Lieb-Robinson bounds consisting of a function which is exponentially growing in time, multiplied by a function 
which decays slower than exponentially in space. As a result, these types of decay do not give rise to a Lieb-Robinson 
velocity; for example, a function exp(t)/Z'^ is only small for t which is logarithmically large in /, while a function 
exp(t) exp(— -v/Z) is only small for t of order Vl. 

A further important thing to remember is that the sum over sites in Eq. (|27p can worsen the bounds beyond what 
one might expect. That is, the decay in the error bound may be not as rapid as the decay in the interactions. For 
example, if we have a two-dimensional lattice, with an interaction H that is a sum of Hz with each Hz having \Z\ — 2 
(i.e., each term in the interaction acts on only two sites), with ^ l/diam(Z), then we find that if the distance 

between a site i and a set Y is equal to I, then the first term in Eq. (|27|) is of order l/l. The second term is of 
order ^^(l/dist(i,j))(l/dist(j, y)) where the sum is over sites j. However, one may see (by replacing the sum by 
an integral) that the second term diverges logarithmically as the lattice size tends to infinity. So, in fact in this case 
no Lieb-Robinson bound can be proven. One needs a faster power-law decay than this to prove the Lieb-Robinson 
bound for a two-dimensional lattice. General consideration of different power-law decays was given in 

One way of treating these was considered in [g^ . We define 

Definition 1. A function K{1) is reproducing for a given lattice A if for any pair of sites i,j we have 



^K{A\si{i,m))K{Aisi{m,j) < XK{dist{i, j)), 

m 



(31) 
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for some constant A. 

For a square lattice in D dimensions and a shortest-path metric, a powerlaw K(l) ^ is reproducing for sufficiently 
large a. An exponential decay is not reproducing. However an exponential multiplying a sufficiently fast decaying 
power is. Using this definition and Eq. (P7)). suppose that \\Hz\\ < K{dia,m{Z)) for some reproducing K. Then, we 
find that Eq. (gT]) is bounded by 

2K{dist{^, Y))i2\t\ + + ^^A^ + ... (32) 

< 2i.(dist(z,r))^^^HM^. 

A 

Note, as mentioned above, that if K decays slower than exponentially, then this function does not lead to a Lieb- 
Robinson velocity; the time t at which we have a meaningful bound will grower slower than exponentially in dist(«, Y). 

A useful trick is that if we start with a function which is not reproducing (such as an exponential) in many cases 
we can bound it by a slightly slower decay exponential times a sufficiently fast power law to arrive at a reproducing 
decay function. 



III. CORRELATION DECAY 



We now consider the decay of correlations in the ground state of a Hamiltonian with a spectral gap and with a 
Lieb-Robinson bound. Such a decay was first proven in [3]. In addition to the interest in this result for itself, we 
consider it because it introduces a combination of two techniques which will be particularly useful. We combine 
Lieb-Robinson bounds with a set of tools using the Fourier transform. Physically, one may imagine that a spectral 
gap AE sets a time scale, AE^^, and then the Lieb-Robinson bound allows us to define a length scale v^jiAE^^ . 
The techniques we introduce combining the Lieb-Robinson bound with the Fourier transform make this statement 
precise. 

The statement of the theorem with a unique ground state 3 is 

Theorem 2. For a quantum lattice system with a unique ground state and a spectral gap AE, and any operators 
Ax, By supported on sets X,Y, we have 

{iPo,AxBYt^o) ~ {^o,Ax^o){-^o,By^o)\ < c{exp(-ZAS/2z;iflJ + min(|A|, |y|)g(0}||OA||||OB||, (33) 

for some constant C, where I = dist(A, Y). We can also consider systems with multiple ground states 4'o as above, 
with a spectral gap AE separating those states from the rest of the spectrum. In this case^], we can instead obtain 
bound the quantity (iPq^AxByiPq) — {ip^, AxPaByipQ, where 

^'o = EW)(V'ol (34) 

a 

is the projector onto the ground state subspace. The bound on this quantity {ipQ, AxByipo) — (?/'" j ^^^o-Byi/jQ is 
equal to the right-hand side of Eq. (j33p plus an additional term which is proportional to the energy difference between 
the ground states, and in fact this result reduces to Eq. (|33p in the case of a single ground state. 

Before sketching the proof of the theorem, we discuss the application of the theorem in various settings. Consider 
the transverse field Ising model Hamiltonian in the paramagnetic phase, B » J, with spin-1/2 on each site. In this 
case, there is a unique ground state and a spectral gap, so that all correlations decay exponentially in space. On 
the other hand, we may consider the Hamiltonian in the ferromagnetic phase, J >> B. For B — 0, the model has 
two exactly degenerate ground states, corresponding to all spins pointing up or all spins pointing down. In fact, the 
Hamiltonian has a particular symmetry. It commutes with the unitary operator Hjl^'^f )■ This operator ffips all the 
spins. So, the eigenstates of the Hamiltonian can be chosen to be eigenstates of this particular unitary. To do this, 
at _B = we may choose the ground states to be symmetric and anti-symmetric combinations of the states with all 
spins up or down: 

^± = ^(itn->± I ;;;..•>). (35) 



For B > but B still much smaller than J, there starts to be a small splitting (exponentially small in system size) 
between the two lowest energy states, and the gap to the rest of the spectrum remains open. Now, consider a local 
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operator such as Sf. We will see what the action of this operator is in the ground state sector. That is, we will 
construct a 2-by-2 matrix M whose matrix elements are the matrix elements of 5*^^ in the ground state. For B = 0, 
we have 



(^+,^fV+) = (V'-,^f^-)=0, (36) 

and 

(V;+,^f^-)-(V-,Sf^+)-l/2, (37) 

so 

(33) 



1/2 



For i? > 0, we instead find that 



where m is the "order parameter" (m is greater than zero in the ferromagnetic phase and m < 1/2 for B > 0). 

The nontrivial matrix elements of in the ground state sector are what give rise to the long range correlations: 
the theorem above shows that {4'+, Sjtfj+) approaches when dist(i,j) gets large. This behavior will contrast 
strongly with the topologically ordered case below. 

We only sketch the proof in the case of a unique ground state (the reader is invited to consider the case of multiple 
ground states). Define By to be the positive energy part of By- That is, let be a basis of eigenstates of H, with 
ipi having energy Ei. In this basis, we define By to have matrix elements 

(S+)., = {By\,,6{E, - E,), (40) 

where Q{x) is the step function: 0{x) = 1 for x > 0, 0(0) — 1/2, and 9{x) = for a; < 0. Without loss of generality, 
assume that {ipQ, Axtpa) — (V'Oj -ByV'o) = 0. Then, we find that 

(V-o, A^SyVo) = {^Po,AxB+i,o) = (^0, [Ax,B+]i,o), (41) 

since ByV'o = Byipo as there are no states with energy less than ipQ. Our strategy is to construct an approximation 
to By, which we call By, which has two properties. First, By has small commutator with Ax- Second, 

B+\^o) ^ B+\i;o), (42) 

and 

{^o\B+ « (^o|S+ = 0, (43) 

where the error in the approximation is discussed below. Then, by making the error in the approximations small, and 
by making the commutator [v4x,i3y] small in operator norm, we will be able to bound ("00; [^x,5y]V'o) as follows: 
we first show (using (|42l43p ) that {tpo, [Ax, iJyJi/'o) is close to {ipo, [Ax, Byjipo). Then, we show that that quantity is 
small using a bound on the operator norm of the commutator. 
To do this, we define By by 

B+^^ lim / dtBy(<)-^exp[-(iAi?)V(2(7)], (44) 

where q is a constant that we choose later. Note that as q — > cx), we find that the Gaussian on the right-hand side of 
Eq. (|44|) gets broader, and By converges to By. So, taking q large will make it easier to satisfy (|42I43|) . Conversely, 
for q large, we will be able to show that || [Ax, By] || is small, since for large times t the integral over t in pi)) is cut off 
by the Gaussian, and for short times t the commutator of Ax with By(t) is bounded by the Lieb- Robinson bound. 
Using the energy gap, one may show that 

B+[i^o) - B+[i;o) <Cexp(-q/2)||Sy||. (45) 
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FIG. 1: Plot of a step function (solid line) and a sketch of a step function convolved with a Gaussian (solid line). Dashed line 
at Ai? shows that the difference between the functions is small at sufficiently large frequency. 



To show this, we bound the absolute value of the matrix element (V'i, (By — By^-\\)Qj for some i > 0. This is equal to 



|(V'»,(B?-B?)V'o)| = |(Sy),o| 1- / dtexp[i(£;,-£;o)t]^— exph(tAS)7(2<z)] 



(46) 



The integral in the above equation is equal to the Fourier transform of the function exp[— (tAi?)^/(2(j)]. This 
Fourier transform is the convolution of the Fourier transform of the step function with a Gaussian. In the figure, we 
sketch the Fourier transform of this function, as well as the step function. For a narrow width of the Gaussian in fre- 
quency (which occurs if g is large), the Fourier transform of the two functions are very close for energy above Ai?. One 
can bound the difference in the Fourier transform by Cexp(— (jf/2). So, |(i?y — i3y-l-)'0oP < X]i>o l(^i')ioP exp(— g), 
so Eq. (j45|) follows. Similarly, one can bound the quantity in (j43p in the same way. 
One may also bound || [Ax, by a triangle inequality as 



< — lim 

27r £^0+ 



dt\\[Ax,BY{t)] 



1 



it 



exp[-{tAE)y{2q)]. 



(47) 



To bound this integral, we split the integral over t into times t less than //wl_r and t > I/vlh. For t < 1/vlr, we use 
the Lieb-Robinson bound to bound the integral by a constant times (7(Z)|X|||Ax||||i?y||. For t > 1/vlr, we use the 
Gaussian to bound it by a constant times 



\Ax\\\\BY\\cM-{lAE/vLR)y2q]. 



(48) 



We now pick q to minimize the sum of terms in (j45l48p . Picking q = IAE/vlh is the best possible choice, and gives 
the error bound in the theorem. 

The most important thing to remember from this sketch is: we combine Lieb-Robinson bounds and Fourier trans- 
forms. We have some function which is not smooth (like the step function). We approximate it by some smooth 
function (the step function convolved with a Gaussian) to obtain a function whose Fourier transform decays rapidly 
in time (in this case as 1/t times a Gaussian). This then allows us to apply the Lieb-Robinson bounds. The same 
ideas are behind the technique of quasi-adiabatic continuation later. 



IV. TOPOLOGICAL ORDER 



There are many different properties that characterize "topological order" . Consider a Hamiltonian such as Kitaev's 
toric code modelfS]. We will not review this model here, since it is well-explained elsewhere. However, this model 
has several unique properties. On a torus, the model has 4 exactly degenerate ground states, with an energy gap to 
the rest of the spectrum. Surprisingly, however, on other topologies, the ground state degeneracy is different. For 
example, the model has a unique ground state on a sphere. This contrasts strongly with the case of the transverse 
field Ising model mentioned above, which has either 1 ground state (in the paramagnetic phase) or 2 ground states 
(in the ferromagnetic phase) independent of the topology of the lattice. This property of the ground state degeneracy 
depending upon topology is one characteristic of a Hamiltonian with topological order. Another characteristic is 
certain corrections to the entanglement entropy d US]. We will not consider these corrections here. 



10 



In fact, we would like to regard both of these properties (the ground state degeneracy and the entropy) as being 
secondary to one particular defining feature of topological order which we now explain. Suppose a model has k 
different ground states, "00 ^'^^ ~ Then, given any operator O, we can project this operator O into the 

ground state sector, defining a fc-by-fc matrix M whose matrix elements are 

M,b = (Vg,OV'g). (49) 

In the case of the toric code on a torus, the matrix M is a 4-by-4 matrix. Surprisingly, if the operator O is local in 
the sense that the diameter of the support of O is sufficiently small compared to L, then the matrix M is equal to 
a constant times the identity matrix. In the case of the toric code, in fact the matrix Mab is equal to a multiple of 
the identity matrix so long as the diameter of the support of O is less than i/2, so one can consider operators O 
which act on a very large number of sites, and yet still the operator acts just as the identity operator when projected 
into the ground state sector. This behavior contrasts strongly with the case of the transverse field Ising model in the 
ferromagnetic phase, for which the operator acts nontrivially on the ground state sector as discussed above. 

If we slightly perturb the toric code, the model remains in the same phase[lll,[l3|- The 4 ground states will no longer 
be exactl y d egenerate, but the difference in energy between them will be exponentially small as a function of system 
size (see |lll. Il2l| for general upper bounds on the energy splitting). Further, this topological order property will be 
slightly weakened: we will instead have a property that if operator O is supported on a set of diameter sufficiently 
small compared to L (for example, a diameter of at most L/2 will suffice for the toric code), then the corresponding 
matrix M is exponentially close to a multiple the identity matrix (again, see [ill IT^ for general upper bounds) . 

This property, that local operators are close to a multiple of the identity matrix when projected into the ground 
state sector was identified in the Hall effect and in other topologically ordered states [ij. Note that this property 
in fact ensures that the splitting between the ground states is small: the Hamiltonian itself is a sum of local operators, 
so by this property of topological order, the Hamiltonian has almost the same expectation value in each of the different 
ground states. 

We can quantify this topological order as in @. We say that a system has (Z, e) topological order if, for any operator 

supported on a set of diameter as most the corresponding matrix M is within e of a multiple of the identity. 
That is, for some complex number z, we have ||M — zl\\ < e. 

One interesting property of this viewpoint about topological order is that it is a property of a set of states, 
rather than a property of a Hamiltonian. Given any set of orthogonal states, -00 j '^^ '^^'^ ^^k whether these states 
are topologically ordered, independent of whether or not the states happen to be the ground states of some local 
Hamiltonian. 

Note that it is not the case that every state V'o has some partner state tAq which gives us a pair of topologically order 
states. For example, consider a system of spin-1/2 spins and consider the state with all spins up, \i>}^) = \ ftt ■••)• 
There is no state |V'o) ttiat is orthogonal to such that the two states V'o for a = 1,2 are topologically ordered. 

A. Topological Order Under Time Evolution 

Note that we can use the Lieb- Robinson bound to describe the behavior of topological order under time evolution's'] . 
Suppose a state </) has (Z,e) topological order for some Z,e. Thus, there exists another state </>' which is the partner 
of state (j). Let us evolve </> for time t under some Hamiltonian H to obtan "0 = e'xjp{—iHt)(j) (one can consider also 
time-dependent Hamiltonians). Define ip' = exp{—iHt)(l)'. We now wish to show thhat ipyip' retain some memory of 
the topological order in (j), (f>': the length scale I will be smaller and the error e will be larger in a way that we can 
bound quantitatively. 

Let O be any local operator supported on a set of diameter / — m, for some m < v^jit. We project this operator O 
into the two dimensional space of states spanned by 0' and The result is equal to the projection of the operator 
exp{—iHt)0 exp{iHt) into the space of states spanned by </>' and (p. However, the operator exp{—iHt)0 cxp{iHt) can 
be approximated, by the Lieb-Robinson bound, by an operator supported on a set of diameter diam(O) + m < I up 
to an error g(m)|X|||^j!(:|| as given by Eq. ([30)) . We bound \X\ by some constant times l"^ for a d-dimensional lattice. 
Thus, we find that ■0, tp' have (Z — m, e + g{m) const. l"^) topological order for any m < VLRt- Thus, topological order 
cannot be completely destroyed in a short-time. For example, suppose that the interactions in the Hamiltonian are 
such that g{l) decays exponentially in I, and suppose that the initial states 0, 0' have (/, e) topological order with 

1 — L/2 and e exponentially small in L. Then, choosing m = i/4, we find that for times up to L/IAvlh), the states 
(j), (j)' have {L/A, e') topological order, with e' still being exponentially small in L. 

Conversely, topological order also cannot be produced in a short time. If we start with a state such as \iPq) — \ ftt •••) 
and evolve for time t under some Hamiltonian H, we cannot produce (Z, e) topological order for I large compared to 
VLRt and e small. 
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V. FOURIER TRANSFORMS 



As seen by our analysis of correlation functions, Fourier transforms play a key role in the application of techniques 
of Lieb-Robinson bounds to quantum many-body systems. In this section, we collect a few useful facts about Fourier 
transforms. 

First, the Fourier transform of a Gaussian is a Gaussian. There is an "uncertainty principle" at work here: the 
narrower the Gaussian in time, the wider the Gaussian in frequency, and vice- versa. 

Second, the Fourier transform of the product of two functions is the convolutions of their Fourier transform. This 
was used,for example, to estimate the term exp(— (7/2) in the calculation of correlation function decay. 

While we have noticed the general principle that there is a tradeoff between the spread of a function in time and in 
frequency, it will be useful in some applications to have functions which have compact support in frequency. That is, 
we would like to find a function g{uj) which vanishes for > 1, with g(0) — 1, and such that the Fourier transform 
g(t) decays as rapidly in time as possibly. In the classic paper [15], it is shown how to construct such functions g{t) 
such that 



\g{t)\ < OieM-m\m, (50) 
for any monotonically decreasing positive function e{y) such that 

^d, (51) 

y 

is convergent. Further, it was shown that this is the optimal possible decay. For example, the function e{y) may be 
chosen to be 

e(y) = l/log(2 + 2/)2. (52) 

Thus, this function g{t) has so-called "subexponential decay" |T6j. A function f{t) is defined to have subexponential 
decay if, for any a < 1, \ f{t)\ < Ca exp(— t"), for some Ca which depends on a. 

Thus, while we cannot quite obtain exponential decay in time and compact support in frequency, we can come very 
close to it, obtaining functions which "almost decay exponentially". 

We will find it useful also in the next section to have a function F{uj) such that F{iL>) = l/w for > 1 and such 
that F{Ld) is odd and the Fourier transform F{t) decays rapidly in time. We can do this using the functions g{t) 
above as follows. First, we assume without loss of generality that the function g{t) above is an even function of t (if 
not, simply take the even part of the function). Then, define the even function f{t) by 

fit) ^ Sit) ~ git), (53) 

where Sit) is the Dirac (5- function (note that fit) is thus a distribution rather than a function). Thus that the Fourier 
transform /(w) has the property that /(w) = for a; = and /(w) = 1 for > 1. Then, define Fit) by 

Fit) ^'-J dufiu)signit - u), (54) 

where sign(f — u) is the sign function: sign(i — u) = 1 for t > u, sign(t — u) = — 1 for t < u, and sign(O) — (since we 
convolve /(w) against sign(t — u), the resulting F(t) is a function, rather than a distribution). We now show the time 
decay of Fit) and we show that the Fourier transform F(aj) is equal to — 1/w for |a;| > 1, as desired (this calculation 
is directly from [6|). 

Lemma 1. Let Fit) be as defined in \54\ Let -F(w) be the Fourier transform of Fit). Then, 

/•oo 

\Fit)\ < I / fiu)du\, (55) 



Fiuj) = — />). (56) 



12 



Proof. Assume, without loss of generality, that t > 0. Then, we have \F{t)\ < \ /(u)du|/2+ | /(u)du|/2. Since 
/(O) = 0, we have | f{u)du\ = \ /,°° f{u)du\. Thus, \F{t)\ < \ /,°° f{u)du\. 



We have 



Integrating by parts in t, we have 



F{uj) = ^ I dtexp{iujt) I duf{u)sign{t-u). (57) 



F{uj) = dtexp{iujt) J duf{u)S{t-u) (58) 

Note that linit_j±oo(^ / duf{u)sign{t — u)j = 0, so the contributions to the integration by parts from the upper and 
lower limits of integration vanish. □ 

Finally, given the decay of f{t), it follows that F{t) decays subexponentially also. 

A side note: we can also use this idea of compact supported functions in the correlation decay calcuation done 
previously. Using the functions g{Lo) described above, we can construct, for example, a family of functions /(w, e) 
such that lim^^Q f (uj , e) = 1 for lo > 1 and J{uj,e) = for uj < —1, and with the Fourier transform g{t,e) decaying 
subexponentially at large times. We construct this family by taking the Fourier transform of the function g{t)/{it + e); 
the Fourier transform of this function converges, as e — ^ 0, to the convolution of the Fourier transform of f{t) with 
a step function. This family of functions still is singular in the limit e — (a singularity like l/{it + e), just as 
we encountered in the Gaussian function exp[—{tAE)'^/2q\/{it + e) in the correlation decay calculation). Such a 
singularity is in fact unavoidable given the large w behavior of the function. This approach will not give quite as tight 
bounds on the correlation decay (the bounds will be subexponential rather than exponential), but the calculation is 
a little simpler since we will have only error terms involving the bound on the commutator || [Ax, Sy]]!, while with 
these compactly supported functions the difference Byi^o — Byipo will vanish. 



VI. QUASI-ADIABATIC CONTINUATION 



We now consider the problem of how the ground state of a local Hamiltonian changes as a parameter in the 
Hamiltonian is changed. Suppose we have a parameter dependent Hamiltonian, iJ^, where s is some real number. 
Suppose that Hg = J2z Hz(s), with Hz(s) being differentiable. Suppose further that we have uniform bounds on the 
locality properties of the Hamiltonian (for example, for all s, the exponential decay ([9]) holds or some other similar 
assumption holds uniformly in s). The main idea of this section can be summarized in a single sentence as follows: if 
such a Hamiltonian Hs has a lower bound on the spectral gap which is uniform in s and has a unique ground state, 
"00 (s), then we can define a Hermitian operator, called the quasi-adiabatic continuation operator which is local 
(in some slightly weaker sense), such that either dgipois) — iVsipois) or dsijjo{s) ~ iT^si^ais) (whether we want exact 
or approximate equality depends upon the application). We will also present a generalization to the case of multiple 
ground states, and use this result to prove a Goldstone theorem. Finally, we will use these ideas to discuss what we 
mean by a "phase" of a quantum system, and to explore the stability of topological order under perturbations. 

In the next section, we will the quasi-adiabatic continuation operator defined in this section to evolve states 
along paths in parameter space such that the Hamiltonian does not necessarily have a spectral gap for s > 0. This 
continuation along these paths will allow us to prove a higher dimensional Lieb-Schultz-Mattis theorem. The fact 
that we continue along paths which might not have a gap is in fact essential to the proof of that theorem: the results 
in this section show that continuing along a gapped path implies that one remains in the ground state, while our goal 
in the next section is to construct a state which is different from the ground state but still low energy to prove a 
variational result. 

We now define the quasi-adiabatic continuation operator: 

Definition 2. Given a parameter- dependent Hamiltonian, Hs, an operator O, and function F{t), we define the 
quasi-adiabatic continuation operator to he the operator 'D{Hs,0) defined by 

iV{Ho,0) = J F{AEt)exp{iHst)Oexp{-iHst)dt, (59) 

where F is an odd function of time so that V is Hermitian. 
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Note that since F is odd, its Fourier transform F{uj) obeys 

F(0) = 0, (60) 

which will be useful in discussions of Berry phase later. 

Given a parameter dependent Hamiltonian Hg = ^zis), we define 

V,=V{Hs,dsHs). (61) 

We also sometimes write — 'D{Hg,dsHz{s)), so that 

2?,-5]pf. (62) 
z 

We will use two different types of functions F{t) in the definition of the quasi-adiabatic continuation operator. The 
first type of function is the function F{t) constructed in the previous section, such that the Fourier transform of F{t) 
obeys F{uj) = —1/uj for |a;| > 1. This will give, as we now show, 

dsipo{s) = iVsikois). (63) 

We call the quasi-adiabatic continuation operator arising from such a function Fit) an "exact quasi-adiabatic contin- 
uation operator". The second type of function F{t) will lead to only approximate equality 

dsMs) ~ iVsMs), (64) 

and we call this the "Gaussian quasi-adiabatic continuation operator" . 
Let us first show 

dsipo{s) = i'Ds^Jois) (65) 
in the case of an exact quasi-adiabatic continuation operator. We have 

iVsMs) = J F{AEt)exp{tH,t)(d,H,^ eM-tHst)dtMs) (66) 
= I]lV'.(s))(^»(s)l / FiAEt)eM^Hst)(^dsHs)cxpi-iHst)dtMs) 
= J2\M^))iMs)^{9sHs)Ms)) f F{AEt)eMKE^is) - Eo{s))t]dt 



E^;(;y^l^.(.))(^.(.),(a.77.)^o(.)) 

9s'0o(s), 



where tpi{s) for i > denote excited states of the Hamiltonian with energy Ei{s). The second line follows by inserting 
the identity as |V'i(s))('0i('S)|, and noting that propertv [501 implies that the term with i = is absent. The third 
line follows by using the fact that the i^sis) are eigenstates. The fourth line follows from the properties of F{t). The 
fifth line is ordinary perturbation theory. 

In the case of multiple ground states, one can generalize this result as follows. We instead find that if we have 
ground states V'ol*) which are not exactly degenerate, then 

a,Vo (s) = ^T^sMs) + Qabi^li-'^). (67) 

where Qab are the matrix elements of some anti-Hermitian matrix Q (if the ground states are exactly degenerate, then 
ds'4>Q{s) may be ill-defined). There is an important Berry phase property: if the ground states are exactly degenerate, 
then the Berry phase arising from the quasi-adiabatic evolution is the same as the usual non-Abelian Berry phase fisj. 
while small corrections to this result occur if there is ground state splitting. Similarly, if Po{s) is the projector onto 
the ground state sector of Hg, then 



dsPo{s)^i[V,,Po{s)]. 



(68) 
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The other important property that the quasi-adiabatic continuation operator, X's has, in addition to (j66|) . is that it 
is local. Let us assume that is bounded by a constant times ||iJz(s)|| to fix a normalization on how rapidly 

the Hamiltonian changes. Then if the original Hamiltonian H has a superpolynomial decay in its interactions (so that 
decays superpolynomially in diam(Z)), and the lattice is finite dimensional, then 'Dg also has a superpolynomial 
decay: we can write I?s = X]z^z(s), where Dz{s) is supported on Z, with ||-Dz|| decaying superpolynomially in 
diam(Z). Note that Dz{s) is not the same thing as . Indeed, I?f is not supported on Z. 

To prove the locality of Dz, one uses the Lieb- Robinson bounds and the superpolynomial decay of the function 
F(t). Before sketching the proof, let us give the basic idea. Consider a given . in integral in Eq. (|59l) is small, while 
at short time we can approximate exp(iHst){dsHz(s)) exp(—iHst) by an operator supported near Z. More precisely, 
we will decompose Vf = X^z^o ^li^)^ where Oi{Z) is supported on the set of sites within distance / of Z as follows. 
We define 

Oo{Z) F{AEt)d,Hz{s))°{t)dt, (69) 

where, following Eq. pT|) . {dsHz{s))'^{t) denotes an approximation to {dsHz{s)){t) = exp{iHst){dsHz{s)) exp{~iHst) 
which is localized on set Z. We define, for I > 0, 

OiiZ) = J F{AEt)(^{d,Hzis)yit) - {d,Hzis)y-\t)yt. (70) 

Summing over I recovers the desired result. Now, we define Dz{s) to be the sum over Y of the OiiY) which are 
supported on Y . 

As mentioned, one can also consider Gaussian quasi-adiabatic continuation operators. These were the first type of 
quasi-adiabatic continuation operators considered while the exact operators will considered later|T9(. The Gaussian 
quasi-adiabatic continuation operators in some cases lead to tighter bounds. For example, in the Lieb-Schultz-Mattis 
theorem, they lead to tighter bounds than the exact operators due to the faster time decay. On the other hand, in 
many cases the exact operators are much more convenient. At this point, we make a deliberate choice due to the 
nature of these lecture notes. Rather than consider the Gaussian operators in detail (which leads to an enormous 
number of triangle inequalities in the actual calculations, potentially obscuring the physics), we will only use the exact 
quasi-adiabatic continuation operators. This will lead to slightly less tight results in many cases, but the improvement 
in clarity (and generality in considering topological phases later), seems well worth it. 

A. Lieb-Robinson Bounds for Quasi- Adiabatic Continuation 

One particular advantage of considering the exact quasi-adiabatic continuation operators is that we have a Lieb- 
Robinson bound for them. We have shown that the norm of the terms Dz{s) decays superpolynomially in diam(Z) 
(indeed, it decays subexponentially if the original Hamiltonian is a finite dimensional lattice with exponentially 
decaying interactions). This implies (see the previous discussion on reproducing functions) that we have a Lieb- 
Robinson bound for quasi-adiabatic evolution; that is, if we evolve an operator O under the equation of motion 

dsO{s)=i[Vs,0], (71) 

we can prove a Lieb-Robinson bound for 0{s). This is a particular advantage compared to the Gaussian case, where 
the proof of Lieb-Robinson bounds for quasi-adiabatic evolution is much more difficult and the bounds are weakerflT}. 

Note that, as discussed previously in subsection (III C|) . we do not actually obtain a finite Lieb-Robinson velocity. 
That is, if we want to bound the commutator || [Ax(s), , where X is supported on X and B is supported on Y , 
the largest value of s for which we obtain a meaningful bound grows slower than linearly in dist(X, Y). This is not a 
problem in most applications, since in general we will be considering path lengths of order unity, while we will often 
consider distances between sets X, Y which are of order system size. 

B. Goldstone's Theorem 

We now present an application to a non-relativistic Goldstone theorem. This theorem is perhaps not that surprising, 
but the results here (originally in [isj ) are more general and simpler than previous nonrelativistic results [20j. 

Above, we have described the clustering of correlation functions, proving that the connected correlated function, 
{AB) — {APqB), of two operators A,B with support on sets X,Y is exponentially small in the distance dist(A, F). 
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Here, for a system with k different ground states. (O) = k ^ Sa(V''oi ^V'o)- For ^ parameter dependent Hamiltonian, 
we define (O), = k-^Ea{^ois),04'Sis)). 

The goal now is to prove (or at least sketch the proof of) a stronger statement about the decay of correlation 
functions in gapped systems with a continuous symmetry showing that the expectation value {APqB) is small also. 

Goldstone's theorem is a statement that a system with a spontaneously broken continuous symmetry has gapless 
excitations. We first need to define a continuous symmetry, or equivalently, a conserved charge. This means that 

Definition 3. We say that a lattice Hamiltonian H has a conserved charge if the following holds. For every site i, 
there is an operator qi supported on site i, with qt having integer eigenvalues. Let Q = "Ylii^i- Then, we require that 

[Q,H]=Q. (72) 

Further, we assume that \\qi\\ < qmax for some qmax (this is a technical point, needed in the later bounds). Our 
non-relativistic Goldstone's theorem will be the contrapositive of the usual statement of Goldstone's theorem: we will 
show that the presence of a gap (between a degenerate ground state sector and the rest of the spectrum) bounds the 
correlation functions. 

For any set X, we define R{0,X) = J^^^^ exp[igi0]. We consider operators (^Xt'\>y with support on sets 
X^Y which transforms as vectors as follows under this U{\) symmetry: R{—6,X)(j)xR{0,X) = exp[id](l)x and 
i?(-6i, r)0yi?(6i, r) = exp^i^y. 

For example, in a Bose system with conserved particle, the qi can represent the particle number on a given site and 
the operators 4>x, 4>y oan represent creation and annihilation operators for the bosons. For a spin system, the qi can 
represent the z component of the spin on a site and the (px , 4>y '^''^'^ represent raising and lowering spin operators on 
sites. 

We do not require the states in the ground state sector to be degenerate with each other, simply the existence 
of a gap between that sector and the rest of the spectrum. This result is stronger than that in [3] as it is valid in 
arbitrary dimension; it is also stronger than other previous results [20j which either required a unique ground state or 
else assumed an ergodic property which is equivalent to requiring the vanishing of the matrix elements in the ground 
state sector in which case the decay or correlations becomes equivalent to clustering. 

We will show that, for a local Hamiltonian on a finite dimensional lattice with a gap /S.E between the ground state 
sector and the rest of the spectrum that {4)x4>y) is superpolynomially small in dist(X, Y) (in 18], stronger exponential 
results are obtained). To show this, we define a set of parameter dependent Hamiltonians Hg as follows. Let X' denote 
the set of sites i such that dist(X, i) < dist(X, Y)/2, as shown in the figure. Then define Hg = R{X', 0)HR{X', —9). 
Clearly, then, as R{X' , —6) is a unitary transformation, T-Lg has the same spectrum oi% and the ground states oiT-Lg 
are given by VS(^)) = R{X',9)M0)). 

Thus, 

dg{(j)x4>Y)e (73) 

= dg{R{x' ,-e)(j3x'4>YR{x' ,9)) 

= dge-x:p[i9]{(j)x4>Y) 

= i{(j)X(j)Y)^ 

where we used the fact that X d X' while FflX' = so that [(/)y, i?(X', 9)] = and where we evaluate the derivatives 
at 6* = 0. 
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Using Eq. ([68t . however, 



de{(f)x4>Y)o 

= -^deTr{PQ{e)MY) 



(74) 



However, recall that Vq is a sum of terms VziO), arising from the different terms dthetaH z{9) . We have Hz{0) — 
RiX', e)HzR{X', -9). However, if Z is a subset of X' or if Z is a subset of the complement of X', then Hz{9) = Hz[o). 
To see this, note that if Hz is a subset of the complement of X' , then Hz commutes with R{X',d). If Hz is a 
subset of X', then Hz{0) — R{A,9)HzR{A,—9), where A is the set of all sites. Since Hz{0) commutes with Q, 
R{A,9)HzR{A,—9) = Hz- Thus, the only terms that contribute to Ve are indeed those where Z intersects both 
X' and the complement of X' . However, the corresponding terms 'Dz{9) have small commutator with (j)x4'Y by the 
locality of the quasi-adiabatic evolution operator: we can approximate Vg by an operator localized near the boundary 
of X' , a distance I from sets X and Y and so the commutator [(Px'Pyt'^o)] can be shown to be superpolynomially 
small after summing over Z. 

It is interesting to note that the assumption of a finite dimensional lattice is necessary in this derivation (it comes 
in when we sum over Z , and is needed to bound the sum of terms by the number of terms times a bound on the norm 
of each term). We sketch a system which is not finite dimensional, and show how a Goldstone theorem may fail in 
this case. Consider a random graph with V nodes each having coordination number 3. Consider a set of V spin-1/2 
spins, with Hamiltonian 



where the interaction matrix J,j equals 1 if the nodes i,j are connected by an edge on the graph, and zero otherwise. 
The interaction is ferromagnetic, so pointing all spins up (or in any other direction) gives a ground state. Further, 
the Hamiltonian is local, using a shorte st p ath metric on the graph to define dist(i, j). However, a random graph of 
this form is typically an expander graph |2l| with a gap in the spectrum of the graph Laplacian, so a spin-wave theory 
calculation [22| gives a gap in the magnon spectrum. Thus, this system has a set of degenerate ground states and a 
gap. However, the spin correlations do not decay, as {Si ■ Sj) = 1/4 for all 



The Lieb-Schultz-Mattis theorem, proven in 1961 [23|, is a theorem about the spectrum on one-dimensional quantum 
spin systems with symmetries. We present the theorem in slightly more general form for theories with a conserved 
C/(l) charge, as considered later by Affieck and Lieb[24,]. 

Consider a one-dimensional Hamiltonian H, with finite-range interactions (one can consider also sufficiently rapidly 
decaying interactions; we do not consider this case in order to make the discussion as simple as possible but we 
encourage the reader to work out what kinds of decay would still allow the theorem to be proven). Assume that the 
Hamiltonian is translationally invariant, with periodic boundary conditions. Let T be the translation operator, so 



Theorem 3. Consider a one- dimensional, periodic, translationally invariant Hamiltonian with finite-range interac- 
tions, with N sites, conserved charge Q and ground state t/'o • Define the ground state filling factor p by 



Assume that p is not an integer. Then, either the ground state is degenerate or the gap between the ground state and 
the first excited state is bounded by 




(75) 



VII. LIEB-SCHULTZ-MATTIS IN HIGHER DIMENSIONS 



[T,H]=Q. 
Then, 



p = (^o,gV'o)/iv. 



(76) 



A£' < const./iV, 



(77) 



where the constant depends only on the strength J of the interactions in H and the range of the interactions. 
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A) 




FIG. 3: A) One ground state of Majumdar- Ghosh model. Circles indicate lattice sites. Light line around circles indicate that 
they are in a singlet. B) Another ground state. 



Let us give some examples of the application of this theorem. Consider the one-dimensional Heisenberg model: 

H = Y,S^■S,+„ (78) 

i 

with spin-1/2 on each site. We identify the charge by 

= + 1/2, (79) 

so that Qi has integer eigenvalues. The Heisenberg model has not just the C/(l) invariance, but instead has a full 
SU{2) invariance (invariance under rotation). So, if the ground state is non-degenerate, then the ground state has 
spin 0. Let us indeed assume that it is true that the ground state has spin (one can also prove this by other means), 
so also the ground state has — 0, which corresponds to p = 1/2, hence p is non-integer. Thus, this model meets 
the conditions of the theorem and so must obey the conclusions: there must be a state within energy of order 1/A^ of 
the ground state. In fact, the lowest energy state has energy of order 1/A^ above the ground state, and corresponds 
to a "spinon" excitation. The model has a continuous energy spectrum in the infinite N limit. 
Another model is the Majumdar-Ghosh model[25j: 

i i 

This model also meets the conditions of the theorem (it again has SU (2) invariance and the ground state turns out to 
have total spin 0). So, it must meet the conclusions of the theorem. However, this model meets the conclusions of the 
theorem in a different way. It has two exactly degenerate ground states and then a gap to the rest of the spectrum. 
One ground state has spins 1 and 2 in a singlet, spins 3 and 4 in a singlet, and so on. The other ground state is 
translated by one, so it has spins 2 and 3 in a singlet, and so on, and finally spins N and 1 in a singlet. So, both 
ground states are products of singlets. If the constant 1/2 is changed to some number near 1/2, then there appears an 
exponentially small splitting between the two ground states, and a gap to the rest of the spectrum, which stil meets 
the conclusions of the theorem. 

A useful exercise for the reader is the following: consider the state which is the symmetric combination of the two 
ground states of the Majumdar-Ghosh model mentioned above, and call this state -00- It is an eigenvector of T with 
eigenvalue -fl. Now, construct the state 4'lsm and verify that it is close to (within distance 1/iV) the state which is 
the anti-symmetric combination of the two ground states and that it is an eigenvector of T with eigenvalue —1. 

Finally, one can also consider spin systems with spin 1 per site. These systems do not meet the conditions of the 
theorem, since if they have total S*^ = 0, then they have integer p. Hence, the theorem does not imply the existence 
of a gap for these systems. In fact, such systems may have a unique ground state and a spectral gap, the so-called 
"Haldane gap" [26]. 

We now prove the one dimensional Lieb-Schultz-Mattis theorem. The proof is variational. First, note that given 
that [H, Q] = 0, without loss of generality we can assume that each term in H commutes with Q. To prove this, let 
H^EzHz. Then, 

H = Y,H'z, (81) 
z 

where 

1 /■^'^ 

H'z=— d9e^pii0Q)Hzexp{-i9Q). (82) 
27r Jo 
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However, each term H'^ commutes with Q. So, from now on, without loss of generahty, we assume that every term 
in H commutes with Q. 



Define a state ^'lsm by 



N 

V-LSAf = (j|exp(27ri-^qj)^Vo- (83) 



One may show that 



{iPlsm,H^Jjlsm) - {^0, HiPo) < const./N. (84) 

To show this, we first show for every Z that 

ii^LSM, HziPlsm) - (V'o, -H'zV'o) < const./Af2^ (85) 

as may be proven using the fact that Z has bounded diameter and that Hz commutes with Q (we encourage the 
reader to work through the detailed proof; the bound on the right-hand side will depend on the diameter of Z). We 
then sum Eq. ([55)1 over Z to arrive at Eq. 

So, we have shown that ipLSM is close in energy to the ground state. We now show that ipLSM is orthogonal to the 
ground state, which will complete the variational proof of this one-dimensional Lieb-Schultz-Mattis theorem. Note 
that we can assume that the ground state is an eigenvector of T (otherwise, since T commutes with H, the ground 
state is degenerate), so that 

T'ipQ = zV'o, (86) 

for some complex number z with \z\ — 1. Now consider Tip^sM- We will show that tpLSM is also an eigenvector of T 
but with an eigenvalue different from z, which will imply that tpLSAi is orthogonal to "001 completing the proof. We 
have: 



N 



T^LSM = T{J[eM^^^j^qj))^o (87) 

N 



{T(nexp(27rz^g,))r-i}T0o 

N 

z{T(j[exp{2mj^q,))T-^}^o 

N 

z exp ( 27r j _ 1 ) ^ Vo , 



where we define go = Qn (recall that we have periodic boundary conditions). Thus, shifting the summation variable 
j by 1 and recalling that qj has integer eigenvalues so exp(27rzgAr) = 1, we have 



N 



T^LSM = z(J]exp(27rz^gj))V'o (88) 



= z 



N 

W exp(27^^^g,))exp(2^^Q/iV)V'n 



N 

Z 

z exp(27rip)-0LSM 



( n exp(27ri-^gj)^exp(27rip)V'o 



In the above equation, we have assumed, without loss of generality, that i/'o is an eigenvector of Q (otherwise the 
ground state degeneracy follows automatically since [Q,H] = 0), so exp{27TiQ / N)tpo = ex-p(2Trip)ipQ. However, since 
we assume that p is not an integer, we find that iJjlsm is an eigenvector T with eigenvalue z exp(27rip)V'LSAf which 
differs from z, completing the proof. 
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The above proof technique simply does not work in two dimensions. Suppose we have a two-dimensional system, 
with size L in each direction which is periodic in one direction. The energy of the state ipLSM is a constant amount 
above the ground state, and does not go to zero as L goes to infinity. The reason is that the energy per site is of 
order l/L"^, while there are a total of sites. In contrast, in the one-dimensional case, there are only L sites. This 
different scaling between one and two dimensions should be very familiar from statistical mechanics: two dimensions 
is the "lower critical dimension" to break a continuous symmetry. 

There is another physical reason why the one dimensional proof fails in two dimensions. In one dimension, there 



are only two possibilities |24|. Either, the system has a continuous spectrum, or, if it has a degenerate ground state 
and a spectral gap (as in the case of the Majumdar-Ghosh model), there is a discrete symmetry breaking, with a 
local order parameter. For example, in the Majumdar-Ghosh model, for any i, the operator Si ■ Si+i has non-trivial 
action in the ground state subspace. In contrast, in two dimensional system, there might also be topological order. 
We might have a system with degenerate ground states but for which no local operator has non-trivial action in the 
ground state subspace (i.e., every local operator is close to a multiple of the identity when projected into the ground 
state subspace). 

Topological order can arise in spin systems; one set of proposals involves the idea that a so-called "short-range 
resonating valence bond" state describes the ground state [2^. Such states are a liquid-like superposition of various 
singlet configurations, that are in many ways very physically similar to dimer models or to the toric code. So, already 
on physical grounds we expect that we will need some new technique to prove a Lieb-Schultz-Mattis theorem in higher 
dimensions. Conversely, we get a nice payoff from such a theorem in more than one-dimension: it will rule out the 
possibility (for certain systems which obey the conditions of the theorem) of having a unique ground state and a 
spectral gap. If we show then that some system does have a spectral gap then either there is ordinary order (some 
local operator has nontrivial action in the ground state subspace) or there is topological order. So, such a theorem 
can be a route to proving topological order. 

We now sketch the higher dimensional proof 0, [13]. This is intended only to be a sketch. The statement is that: 

Theorem 4. Consider a Hamiltonian H , defined on a finite- dimensional lattice with finite interaction range R and a 
bound on interaction strength J. Let H have translation invariance in one direction with periodic boundary conditions 
and have a length L in that direction. Let the lattice have a total of N sites, with N bounded by a constant times a 
polynomial in L. Let H have conserved charge Q and ground state ipQ. Define the ground state filling factor p by 

P= (V'o,Wo)/i. (89) 

Assume that p is not an integer. Then, either the ground state is degenerate or the gap between the ground state and 
the first excited state is bounded by 

AE < const, log(i) /L, (90) 

where the constant depends only on R,J , Qmax, o,nd the lattice geometry. 

The theorem can be extended to sufficiently rapidly decaying interactions also. The bound that N is at most a 
polynomial times L implies that the theorem works for aspect ratios of order unity (for example, an L-hy-L square 
lattice in two dimensions) or even aspect ratios which are quite far from unity (an L-hy-L^ square lattice, for example). 
Note that the bound is slightly weaker than in one dimension (we have log(L)/L instead of l/L). 

Finally, the fact that p = {ipo, Qipo) / L is a minor annoyance in the statement of the theorem. In a spin-1/2 system 
on a square latttice with odd width (i.e., an L-by-M lattice with translational invariance in the first direction and 
with M even) we do indeed find that p is non-integer. However, for such a system on a lattice of even width, the 
theorem does not work. In fact, there are counterexamples to a conjectured theorem with even width (a spin ladder, 
consisting of an i-by-2 system of spin-1/2 spins with Heisenberg Si ■ Sj interactions between nearest neighbor spins 
has a unique ground state and a gap). However, we expect that if a two-dimensional system on an L-hy-M lattice 
has translation invariance and periodic boundary conditions in both directions, then the goes to zero as both L and 
M get large. This has not been proven yet. Still, the theorem above covers a wide variety of cases with a minimal 
number of assumptions (only one direction of translation invariance required). 

Further, the translation invariance in at least one direction is a necessary condition. The reader is invited to work 
out a counter-example if no translation invariance is assumed. 

The sketch of the proof is as follows. It again is variational. It is also a proof by contradiction. That is, we assume 
that the Hamiltonian has a spectral gap AE and use this assumed spectral gap to construct a variational state which 
has low energy. If the initial gap is large enough (larger than a constant times \og{L)/L) the variational state will 
have energy less than \og{L)/L proving the theorem by contradiction. 
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We begin by defining a parameter-dependent family of Hamiltonians, Hg. These Hamiltonians arc defined by 
"twisting the boimdary conditions" in one particular direction, the direction in which the lattice is translation invari- 
ant. Let us label the coordinate of a site i in this direction by x{i), with < x{i) < L. To define the flux insertion 
operator, we need to define the Hamiltonian with twisted boundary conditions. Let Qx be defined by 

l<x{i)<L/2 

Qx= (91) 

i 

where x{i) is the x-coordinate of site i. That is, Qx is the total charge in the half of the system to the left of the 
vertical line with x = L/2 + 1 and to the right of a; = 0. Let 

H{eue2) = Y,Hz{9i,e2), (92) 

z 

where i?z(^i, ^2) is defined as follows. If the set Z is within distance R of the vertical line a; = 0, then Hz{0i,92) = 
exp{i6iQx)Hz exp{—i6iQx)', if the set Z is within distance R of the vertical line x = L/2, then Hz{0i,62) = 
exp(— i^2<3x)^fz exp(i^2<3x); otherwise, Hz{6i,02) = Hz- Note that, 

H{0, -6) = exp{ieQx)H exp{-ieQx)- (93) 

This unitary equivalences implies that H{9, —6) has the same spectrum as H which will be useful below. 

The introdiiction of the two different vertical lines is an important technical trick. We now define an operator Wi 
which generates the quasi-adiabatic evolution along the path where 9i evolves from to 27r and 62 = 0. That is, 
define Vg to generate the quasi-adiabatic evolution for Hs = H{s, 0) and let 

Wi = exp( / deVl), (94) 
Jo 

where the exponential is 9 ordered. Similarly, let W2 generate quasi-adiabatic evolution along the path where O2 
evolves from to 27r and 9i = and let generate the quasi-adiabatic evolution for Hg = H{0, —s). Finally, let W 
generate quasi-adiabatic evolution along the path 9\ = —92 = 9 as 6 evolves from to 27r. 

An important point: we do not assume that the gap remains open along the paths above used to define Wi,W2. 
We simply use the assumed initial gap at 0i = ^2 = and then evolve quasi-adiabatically as if the gap remained open 
along the path. 

Using locality of the quasi-adiabatic evolution operators, if the gap is sufficiently large, one may show that 

W"iW2 « W2W"i « W, (95) 

where the approximation means that W1W2 — W is small in operator norm. The size that the gap needs to be and 
the magnitude of the error ||H^iW2 — W\\ both depend on the quasi-adiabatic evolution operator we use. Roughly, 
the quasi-adiabatic evolution operator is supported near the line x{i) = 0, up to a length scale which is inversely 
proportional to the gap. We need this length scale to be small compared to L/A so that we can approximate the 
operator by an operator supported within distance iy/4 of the line x{i) = 0; we make the same approximation 
for so that in this case, the operators T)^ and can be approximated by operators supported on disjoint sets. 
Using exact quasi-adiabatic evolution operators, one finds that the gap needs to be at least f{l)/L for some function / 
growing slower than any polynomial, while for Gaussian operators one can choose f{L) to be a constant times log(L). 
We consider the variational state: 

WiVo- (96) 

We wish to show that this state has low energy. Note that W\ is unitary. The Hamiltonian H is & sum of terms Hz- 
We will show that each term Hz has roughly the same expectation value in the state Wiipf) as it does in the ground 
state. If Z is far from the line x{i) = (for example, far can mean that the distance is at least L/4), then this follows 
directly from the locality of the quasi-adiatic evolution operator: Hz almost commutes with W\ : 

(W^i^cJ^zWiVo) ~ (V'c^zV'o). (97) 
Conversely, if Hz is near the Inie x{i) — (say, the distance is less than i/4), Hz almost commutes with W2- So: 

{W^^o,HzW^i)o) « {^Ja,WlwlHzW2Wi^o) (98) 
« {iPq,W^HzW^o)- 
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x=0 x=L 

FIG. 4; Twist in boundary conditions are applied at two places, at x{i) = = L (along the boundary of the systems) and at 
x{i) — L/2 (along the dashed line). We want to approximate Wi by an operator supported on the upward slanting grey lines 
(near x{i) = 0) and to approximate the operator W2 by an operator supported on the downward slanting grey lines, so that 
Wi, W2 will approximately commute. This requires a gap sufficiently large compared to 1/L. 



Note, however, that W describes evolution along a path of Hamiltonians which are all unitarily equivalent to H{0, 0) 
so all Hamiltonians along this path have the same spectral gap. Further, H{2tt, — 27r) = -ff (0, 0) So, WipQ = wtpo for 
some complex number w with \w\ = 1 (we will worry about this phase w in the next paragraph, for now it doesn't 
matter). So, {Wi^po, HzWitpo) « (i^o^HzM for aU Z. 

This proves that Wiipo is a low energy state. We now need to show that Wi4>o is orthogonal to the ground 
state. We do this by showing that it has a different expectation value for the translation operator, T, than the 
ground state does. Suppose {ipOjTipo) = z for some complex number z. Then, consider {tl^Q^WlTWiTpo) . This equals 
{'ipo,Wl{TWiT-^)Ti)o) = z{tpo,Wl{TWiT-^)tl)o). This is approximately equal to z{'iPo,wlwl{TWiT-^)W2i^o) ■ 
Note that W^Wl is close to W\ so W1W21P0 is close to w4>o where w is some phase as mentioned above. Similarly, 
W2WI {TWiT^^)W2tpo is close to w'tpo for some other complex number w' with \w'\ = 1. To see this, note that 
(TWiT~^) describes quasi-adiabatic evolution where we twist the boundary conditions along the line x{i) — 1 rather 
than along x{i) = 0, so {TWiT~^)W2 is close to some operator W' which describes quasi-adiabatic evolution of a 
Hamiltonian with twisted boundary condition by 9 along line x{i) = 1 and by —6 along x{i) = L/2. So, {Witpo, TWitpo) 
is close to ww'z. We now just need to work out the phases w, w'. 

However, using the property (1501) . one can show that w = exp(27ri('(/'oj QxV'o))- Similarly, the phase w' depends on 
the expectation value of qi summed over 1 < x{i) < L/2. So, 

w'w = exp(27r«(V'o, ^ qtipo))- (99) 

i,x(i) — l 

This expectation value is non-integer by assumption. So, WiijjQ has a different expectation value for T than ipo, so it 
is orthogonal to ipQ. 

Some heuristic comments on why physically our proof of the theorem is necessarily a proof by contradiction. That 
is, why we assumed a gap at the beginning of the proof. We use the idea of twisting the boundary conditions. A state 
such as an anti-ferromagnet will strongly resist this twist in boundary conditions; that is, the ground state energy will 
change by an amount of order unity when we impose this boundary twist in a two-dimensional system (and by an even 
larger amount in higher dimensions). Thus, for a state such as an anti-ferromagnetic (which is gapless), there is no 
reason to expect that the procedure we described of twisting boundary conditions and following the quasi-adiabatic 
evolution of the state along the path will give any useful results. However, suppose we have a system which has a 
gap. By the theorem we have just sketched, such a system cannot have a unique ground state and then a gap to 
the next lowest energy state. So, we instead want to consider a state with a degenerate ground state and a gap. 
Such a system could be a valence bond solid or a resonating valence bond system, among other possiblities. In such 
a system, twisting the boundary conditions does not lead to a large energy cost. Instead, it typically costs only an 
exponentially small amount of energy to twist the boundary conditions. Then, when we start at one ground state 
and quasi-adiabatically evolve it, twisting the boundary angle 6 from to 27r, we transform it a state close to an 
orthogonal ground state. This is analogous to the discussion in one dimension, where starting with one of the ground 
state of the Majumdar-Ghosh model and constructing the state tpLSM, gave us something close to the other ground 
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state. So, the physical idea of the higher-dimensional proof is that if there is a sufRciently large gap from the ground 
state sector to the rest of the spectrum, then we can use the quasi-adiabatic continuation to construct a unitary that 
transforms one ground state into another. 

VIII. WHAT IS A PHASE? 

What is a phase of a quantum many-body system? We are used to the idea that physical systems have distinct 
phases, with phase transitions between them. For example, water can appear as ice, water, or steam (and further, 
there are many distinct phase of ice). This discussion of the properties of water is a discussion of systems at non- 
zero temperature, while our focus in this notes is on quantum systems at zero temperature, but many of the same 
phenomena occur in both cases. For example, we do not actually consider steam and water to be distinct phases of 
matter. While usually water turns into steam by being boiled (a phase transition, where the energy is non-analytic in 
the thermodynamic limit), we can also move from water to steam without any phase transition, by following a path 
in the two-dimensional plane of temperature and pressure. A similar phenomenon occurs in the transverse field Ising 
model. Suppose we consider the model with an additional parallel magnetic field so that the Hamiltonian is 

N-l N N 

h = -jJ2s^s^+, + bJ2s! + hY,s^- (100) 

i—l i—1 i—1 

Suppose J >> B and H > 0. Then, the system has a unique ground state; at i? = 0, this ground state is the state 
with all spins up, while for B > 0, there are quantum fiuctuations about this state. 

When H changes sign, the ground state changes from all spins up to all spins down in the case B — 0, crossing a 
phase transition. This is a zero temperature phase transition. In the case of i? = 0, this phase transition is a level 
crossing: at _ff = 0, there are two exactly degenerate ground states. For B ^ 0, this level crossing becomes an avoided 
crossing which we now describe. Recall that we said that for S 7^ but H = 0, the system has two ground states with 
an exponentially small splitting between them. The behavior of the two lowest energy states as a function of H can 
be roughly understood in the following toy model. Suppose that H is very small. Then we can focus on just the two 
lowest energy states and use Eq. pQ]) to arrive at the following two-by-two Hamiltonian (we obtain this Hamiltonian 
by projecting the term X]i=i i'^to the ground state subspace): 

where t is some exponentially small splitting between the two lowest states and = |A| is the size of the system. 

By a change of basis (going to symmetric and anti-symmetric combinations of the two ground states), we arrive 
instead at 

The two different basis vectors here correspond to the spin up and spin down ground states. 

While this Hamiltonian is valid for small H, for larger H we need to worry about the excited states. However, 
this Hamiltonian already reveals the essential point, namely that since t is exponentially small as a function of A^, 
in the limit of A'' — >■ 00, the ground states energy per site is a non-analytic function of H. However, this behavior 
is a lot like changing from water to steam by boiling: we can also move from H < to H > Q without crossing a 
phase transition. Instead, one should follow the path of first making H large and negative, then decreasing J until 
J << B, then changing the sign of H (which does not involve a phase transition since J « B) and then increasing 
J, as shown in Fig. (|VIIip . 

We are motivated by this analysis to adopt the following definition of a quantum phase: two Hamiltonians, Hq and 
Hi describe systems in the same quantum phase if both and Hi have a spectral gap, and one can find a smooth 
path Ha connecting Hq and Hi which keeps the gaps open and keeps the Hamiltonian local. (As a technical point 
for those interested, in some cases it may be more appropriate to consider a "stable limit" when describing equality 
of quantum phases, as in the case of topologically ordered phases of free fermion systems (29|]). 

Now, one may choose to have a more refined notion of a quantum phase, which takes into account symmetries; for 
example, one might wish to insist that there is a path connecting H^ to Hi which respects symmetries such as an 
Ising symmetry. We do not consider this kind of restriction on the definition here. Our interest is instead the case of 
systems which are in distinct phases even without any assumptions on the symmetry. That is, Hamiltonians which 
cannot be connected by such a smooth path. 
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FIG. 5: Energy of lowest two states as a function of parallel magnetic field H. This is a sketch. The crossing is an avoided 
crossing but the splitting between states is exponentially small at -ff = 0. At -ff of order there is another avoided crossing 
as the energy gap becomes of order unity. 



We can use the ideas developed in these notes to prove that certain Hamiltonians cannot be connected by a smooth 
path of local Hamiltonians without closing the gap or without the path length being long (where "long" means that 
the path length diverges as the system size gets large). Consider, for example, a toric code on a torus. This has 
four ground states and a spectral gap. Call this Hamiltonian Ho. Consider instead a Hamiltonian consisting of two 
copies of the transverse field Ising model on a torus with J » B. Call this Hamiltonian Hi (if one wants a better 
statement of TJi, break the square lattice into two different sublattices, and have interactions only between spins on a 
given sublattice, so that way we have two copies of the transverse field Ising model with the same number of degrees 
of freedom as in the toric code system). Note that iJo, Hi both have 4 ground states and a gap. So, can we find a 
path connecting Hq to Hil The answer is no. If such a path existed, then we could use quasi-adiabatic continuation 
to evolve the four ground states of Hq to produce some linear combination of the four ground state of Hi. However, 
this would imply that (recall the discussion in subsection IIV Al and the fact that we have a Lieb- Robinson bound for 
the quasi-adiabatic continuation operators) the ground states of Hi are also topologically ordered. Since this is not 
true, no such path can exist. The reason we need to assume that the path is not "long" is that the length of the path 
plays the role of time in the Lieb-Robinson bound, and recall that in the discussion of the behavior of topological 
order under time evolution we only showed that topological order could not appear after evolution for a short time, 
but not for arbitrary time. 

We can use similar arguments to show that Hamiltonians cannot be connected by such a smooth path even when the 
ground state is unique (again, without closing the gap and without the path being long) . One of the key properties of 
systems like the toric code is the particular set of expectation values they have for certain operators they have called 
string operators; these are operators which are products of single-site operators around a loop which are analogous 
to Wilson loops in gauge theories. One can define "dressed operators" by quasi-adiabatically continuing these string 
operators along a path. Thus, given a string operator O for a Hamiltonian Hq, we can define an operator 



H 



J/B 



FIG. 6: Path to follow. Thickened line on axis denotes B/ J less than the critical value. 
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where 

U = S exp(i / dsVs), (104) 
Jo 

where the calhgraphic s in front of the exponential denotes that it is an s-ordered exponentiaL Then, the operator O 
has the same expectation value in the ground state of Hi as O does in the ground state of Hq . This dressed operator 
O is precisely equal to the operator 0{s = 1) as described by the evolution of Eq. ([7T|) . 

Further, if two operators O, O' anti-commute with each other, then the operators O, O' also anti-commute with 
each other. Now, consider the transverse field Ising model in the phase B » J. This has a unique ground state. 
Similarly, the toric code on a sphere has a unique ground state. However, these two models cannot be connected by 
a continuous path. To show this, note that if they were connected, then we could also connect the toris code to the 
transverse field Ising model in the phase B ^ 0, J = 0. However, the ground state of this model with _B 7^ 0, J = can 
be shown to be inconsistent with the properties of such dressed string operators. Here is a sketch (this sketch depends 
on properties of the toric code which are not discussed in these notes and need to be read elsewhere): Consider an 
electric loop operator, indicated as the solid line in (A) of the figure. Call this operator E. Let the dashed line 
represent a magnetic loop operator, which we call M . These operators commute with each and both have expectation 
value 1 in the toric code ground state. In (B) of the figure, we show that the operator E can be written as a product 
of two different operators, which we call Ei and E2 which act on part of the loop. Note that {Ei,M} = {E2, M} = 0. 
Hence, the expectation value of E1ME2 is equal to minus 1 in the toric code ground state. Now consider the dressed 
operators Ei,E2,M assuming the existence of a path connecting the toric code to the transverse field Ising model. 
Since E1E2 would have expectation value 1 in the transverse field Ising model ground state, V'li the state EiE2ijJi 
must be simply the product state of all spins pointing along the transverse magnetic field. However, using locality of 
the dressed operators, Ei can be approximated by an operator supported near the support of Ei. Hence, the state 
i?2V'i niust have all of the spins which are far from the support of Ei pointing approximately along the transverse 
magnetic field. Thus, the state i?2V'i has all of its spins, except those near the upper and lower ends of the support 
of E2 (those close to the support of both Ei and E2) pointing approximately along the magnetic field. Similarly, the 
operator M can be approximated by an operator supported near the support of M; however, this means that M can 
be approximated by an operator, which we call M', which is not supported near the upper and lower ends of i?2- 
As noted, the spins away from the upper and lower ends of the support of E2 are aligned approximately along the 
magnetic field, and those are the only spins in the support of M' . Since Mipi = V'l, we have M'tpi close to ipi and 
so, using the fact that the spins in the support of M' are almost aligned with the field in the state £'2V'i and in the 
state ipi, we find that M'E2iIji is close to E2ipi- Thus, ME21P1 is close to Eiipi- Hence, we find that the expectation 
value of E1ME2 is close to unity in the transverse field Ising model, while it was close to minus 1 in the toric code, 
giving a contradiction. 

One may also choose to take into account other (anti-unitary) symmetries such as time-reversal symmetry. These 
symmetries are well- understood in the non- interacting case [29], but only a limited understanding has been obtained 
in the interacting casefsol]. In general, even without symmetries the classification of phases of matter of interacting 
systems under the definition above is only in the earliest stages. One interesting case is that lattice Hamiltonians 
are known, the so-called "Levin- Wen models" jsij , which realize certain two-dimensional unitary topological quantum 
field theories (TQFTs), in particular those theories which are quantum doubles. Thus, a classification of interacting 
phases of matter requires a classification of these TQFTs. Some results on classification of TQFTs were obtained in 
[33 • However, it is likely that the full classification of lattice models includes many other phases in addition to those 
described by TQFTs, so this classification is a problem for the future. 
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IX. STABILITY OF TOPOLOGICALLY ORDERED PHASE 



In this section, we briefly mention certain recent results on the stability of quantum phases. The analysis throughout 
these notes has always dealt with systems with a gap. In some cases, we were able to prove either the absence of a 
gap or the degeneracy of a ground states (such as in the Lieb-Schultz-Mattis theorem). However, the reader may be 
wondering: how do we know a system has a gap? Similarly, suppose we have proven that a certain model, such as 
a toric code or Levin- Wen model, is in a topologically nontrivial phase and cannot be connected to a topologically 
trivial phase without closing the gap. However, what happens if we slightly perturb the Hamiltonian? Suppose we 
consider a Hamiltonian 

H = Ho + sV, (105) 

where Hq is some unperturbed Hamiltonian describing a topologically nontrivial phase, s is some real number, and 
y is a perturbation. Does the model remain in the same phase for sufficiently small s? Does the gap remain open? 
The interesting question here is to consider the case in which y is a sum of local terms. Thus, we want 

V = Y.^z, (106) 
z 

where the operator norms decay rapidly as a function of the diameter of the set Z (just as we required a similar 

decay on the norms of the terms ||i?z|| in the Hamiltonian Hq). A very elementary result is that for any given system 
size, there is an sq such that for |s| < sq the gap remains open: simply use the fact that for any given system size, 
the norm of the operator V is finite, and the gap for s > is lower bounded by 

A£;(s) > A£;(0) - 2s|ly||, (107) 

where AE{s) denotes the gap as a function of s and the factor of 2 in front of the second term occurs because the 
ground state energy increases by at most s\\V\\ while the first excited state energy decreases by at most the same 
amount. 

So, we may take sq = Ai?(0)/4||F||. However, such a bound, while elementary, is also fairly useless, since it leads 
to an Sq which tends to zero as the system size tends to infinity. Instead, we want a bound which is a uniform 
function of system size. Such bounds were provided in [ill, We will not review them here, except to note 

that using such bounds one can then prove (using quasi-adiabatic continuation) that many of the properties of the 
topologically ordered system (such as ground state splitting, braiding, fusion rules, etc..) remain the same in this 
phase. An interesting open question is to understand the behavior of topological entanglement entropy as a function of 
perturbation. Perhaps some smoothed definition of topological entanglement entropy exists (smoothing over different 
boundaries?) such that it can also be proven to be invariant under perturbations? These problems, and problems like 
the classification of different phase of lattice quantum systems, are problems for the future. 
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